A rate equation approach to cavity mediated laser cooling 



Tony Blake,* Andreas Kurcz, and Almut Beige 
The School of Physics and Astronomy, University of Leeds, Leeds, LS2 9JT, United Kingdom 

(Dated: February 15, 2012) 

The cooling rate for cavity mediated laser cooling scales as the Lamb-Dicke parameter r) squared. 
A proper analysis of the cooling process hence needs to take terms up to rf' in the system dynamics 
into account. In this paper, we present such an analysis for a standard scenario of cavity mediated 
laser cooling with r] <^ 1. Our results confirm that there are many similarities between ordinary 
and cavity mediated laser cooling. However, for a weakly confined particle inside a strongly coupled 
cavity, which is the most interesting case for the cooling of molecules, numerical results indicate 
that even more detailed calculations are needed to model the cooling process accurately. 
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I. INTRODUCTION 

First indications that cavity mediated laser cooling al- 
lows us to cool particles, like trapped atoms, ions and 
molecules, to much lower temperatures than other cool- 
ing techniques were found in Paris already in 1995 [TJ [2] . 
Systematic experimental studies of cavity mediated laser 
cooling have subsequently been reported by the groups 
of Rempe [MS], Vuletic [THin], and others [Ulllj. Re- 
cent atom-cavity experiments access an even wider range 
of experimental parameters by replacing conventional 
high-finesse cavities [13l [14] through optical ring cavities 
[T51 [TB] and tapered nanofibers [T71 [TB] and by combin- 
ing optical cavities with atom chip technology [THl [10] j 
atomic conveyer belts [211 122] , and ion traps f23| . More- 
over, Wickenbrock et al. [Ml recently reported the obser- 
vation of collective effects in the interaction of cold atoms 
with a lossy optical cavity. Motivated by these develop- 
ments, this paper aims at increasing our understanding 
of cavity mediated laser cooling. 

Cavity-mediated laser cooling of free particles was first 
discussed in Refs. [251 [55]. Later, Ritsch and collabora- 
tors |?7ti5T] . Vuletic et al. [31113], and others (51H57] de- 
veloped semiclassical theories to model cavity mediated 
cooling processes very efficiently. The analysis of cavity 
mediated laser cooling based on a master equation ap- 
proach has been pioneered by Cirac et al. [38^ in 1993. 
Subsequently, this approach has been used by many au- 
thors [SnmSj, since the precision of its calculations is 
easier to control than the precision of semiclassical calcu- 
lations. Moreover, cavity mediated cooling is especially 
then of practical interest when it allows to cool particles 
to very low temperatures, where quantum effects domi- 
nate the time evolution of the system and semiclassical 
models no longer apply [28] . 

In close analogy to ordinary laser cooling [ISj-jSn], the 
effective cooling rate of cavity mediated laser cooling 
scales as the Lamb-Dicke parameter 77 squared. A proper 
analysis of the cooling process hence needs to take terms 
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FIG. 1: Experimental setup of externally trapped particles 
inside an optical cavity with coupling constant g and spon- 
taneous decay rates k and F. The motion of the particles 
orthogonal to the cavity axis is confined by a harmonic trap- 
ping potential with phonon frequency 1/. Moreover, a cooling 
laser with Rabi frequency Q is applied. 



of the order 77^ in the system dynamics into account. A 
study of cavity mediated laser cooling which is analogous 
to the calculations presented in Refs. [48lf50] cannot yet 
be found in the literature. It is the purpose of this pa- 
per to present such an analysis for a certain standard 
scenario (cf. Figs. [Tjand 2]). In the following we use the 
master equation approach of Refs. [38], avoid certain ap- 
proximations made in earlier calculations [3111301113], and 
focus our attention on a wide range of experimental pa- 
rameters. For example, our analysis includes the weak 
confinement regime which is of practical interest for the 
cooling of molecules [33H3Z] • 

The reason that our calculations are relatively straight- 
forward is that we replace the phonon and the cavity 
photon annihilation operators b and c by two commut- 
ing annihilation operators x and y. These commute with 
each other and describe bosonic particles that are nei- 
ther phonons nor cavity photons. The corresponding 
system Hamiltonian no longer contains any displacement 
operators and provides a more natural description of the 
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FIG. 2: Level configuration showing the ground jO) and the 
excited state |1) of the trapped particle. Here ojl, i^c, and lvq 
are the frequency of the cooling laser, of the cavity field, and 
of the 0-1 transition of the particle, while 5 and A denote 
detunings. 



cavity-phonon system. Non-linear effects in the interac- 
tion between the vibrational states of the trapped particle 
and the field inside the optical cavity are now taken into 
account by a non-linear term of the form x{y — y^). Us- 
ing the master equation corresponding to this Hamilto- 
nian, we then derive a closed set of 25 cooling equations. 
These are linear differential equations which describe the 
time evolution of x, y, and mixed operator expectation 
values and can be solved analytically as well as numer- 
ically. Applying the same methodology to laser cooling 
of a single trapped particle [50j , we recently obtained re- 
sults which are consistent with the known results of other 
authors piBS] . 

The analytical calculations in this paper confirm that 
there are many similarities between ordinary and cavity 
mediated laser cooling |43j . In the strong confinement 
regime, where the phonon frequency v is much larger 
than the spontaneous cavity decay rate k, we find that 
the optimal laser detuning S^g is close to v (sideband 
cooling). Different from this, one should choose Scs close 
to in the weak confinement regime with <C k in 
order to minimise the final temperature of the trapped 
particle. What limits the final phonon number are the 
counter-rotating terms in the particle-phonon interaction 
Hamiltonian Hi which can increase the energy of an open 
quantum system very rapidly [HIl 122] . 

In addition, we present detailed numerical studies. 
These take all available terms up to order 77^ in the rate 
equations into account, even the ones that appeared to 
be negligible during analytical calculations. Our numeri- 
cal results are in general in very good agreement with our 
analytical results. However, for relatively small phonon 
frequencies v and relatively large effective cavity coupling 
constants Qcb, the stationary state phonon numbers pre- 
dicted by both calculations can differ by more than one 
order of magnitude. This indicates that the analysis of 
the cooling process of a weakly confined particle inside a 
strongly coupled cavity needs to be more precise. Terms 
of order 77'^ should be taken into account systematically 



when modelling the system dynamics with rate equa- 
tions. Notice that the parameter regime where gcs/i^ 
is relatively large is of special interest for the cooling of 
molecules. These usually have a relatively large atomic 
dipole moment but nevertheless cannot be trapped as 
easily as other particles. 

There are six sections in this paper. Section [ll] intro- 
duces the master equation for the atom-cavity-phonon 
system shown in Fig. [T] and simplifies it via an adiabatic 
elimination of the excited electronic states of the trapped 
particle. Section |III| uses this master equation to derive 
a closed set of 25 cooling equations. These simplify re- 
spectively to a set of five effective cooling equations in 
the weak confinement regime and to a single effective 
cooling equation in the strong confinement regime. In 
Section |IV| we show that the phonon coherences and the 
mean phonon number m always reach their stationary 
state. Section fV] analyses the cooling process in more de- 
tail and calculates effective cooling rates and stationary 
state phonon numbers. Finally, we summarise our find- 
ings in Section |Vlj Mathematical details are confined to 
Apps.|XHD] 



II. THEORETICAL MODEL 

The experimental setup considered in this paper is 
shown in Fig. [l] It contains a strongly confined particle 
inside an optical cavity. The aim of the cooling process 
is to minimise the number of phonons in the quantised 
motion of this particle in the direction of a cooling laser 
which enters the setup from the side. Cooling the motion 
of the particle in more than one direction would require 
additional cooling lasers. 



A. The Hamiltonian 

The Hamiltonian of the atom-cavity-phonon system in 
Fig. [1] is of the general form 

H = iJpar + ^^phn + -^cav + -f^L + ^par-cav • (1) 

The first three terms are the free energy of the trapped 
particle, its quantised vibrational mode, and the quan- 
tised cavity field. Suppose, the particle is a two-level 
system with ground state |0) and excited state |1) and 
the energies huo, hv, and hujc are the energy of a single 
atomic excitation, a single phonon, and a single cavity 
photon, respectively, as illustrated in Fig. [2j Then 



-^par 


= hijjQ a a 






-^cav 


= huJc C , 



where the operators a = |0)(1| and cr"'" = |1)(0| are the 
atomic lowering and raising operator, b is the phonon 
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annihilation operator, and c is the cavity photon annihi- 
lation operator with the bosonic commutator relation 



[b,b^] = [c,ct] = 1, 



(3) 



Let us now have a closer look at the two remaining terms 
and i?par-cav in Eq. Q. 
The role of the cooling laser is to establish a coupling 
between the electronic states |0) and |1) of the trapped 
particle and its quantised motion. Its Hamiltonian in the 
dipole approximation equals 



iJL = eD-EL(x,i), 



(4) 



where e is the charge of a single electron, D is the dipole 
moment of the particle, 



D = D 



H.c. , 



(5) 



and El(x, t) denotes the electric field of the laser at po- 
sition X at time t. Moreover, we have 



EL(x,<) = Eoe'(''^-"-"^*)+c.c. 



(6) 



with Eq, Icl, and wl denoting amplitude, wave vector, 
and frequency of the cooling laser. 

The interaction Hamiltonian describing the coupling 
between the particle and the cavity in the dipole approx- 
imation is given by 



-^pai— cav 



eD • Ecav(x) ; 



(7) 



where Ecav(x) is the observable for the quantised electric 
field inside the resonator at the position of the particle. 
Denoting the corresponding coupling constant as g, the 
above Hamiltonian becomes 



par— cav 



= hg{a +cr+)c + H.c. 



(8) 



This Hamiltonian describes the possible exchange of en- 
ergy between atomic states and the cavity. 



B. Displacement operator 

The relevant vibrational mode of the trapped parti- 
cle is its center of mass motion in the laser direction. 
Considering this motion as quantised with the phonon 
annihilation operator b from above yields 



kL • X 



f]{b + b^), 



(9) 



where the Lamb-Dicke parameter i] is a measure for the 
steepness of the effective trapping potential seen by the 
particle W\. Substituting Eqs. into Eq. Q, we 

find that the laser Hamiltonian is a function of the par- 
ticle displacement operator |53j 



D{ir,) 



This operator is a unitary operator. 



(10) 



(11) 



with 



D{ir])bD{iT]y = b + i-q, 
D{i-q)HD(iri) = b-i-q. 



(12) 



Using this operator, Hi^ can be written as 

i/L = e [Doi (J- + H.c] • E; D{iri) e'"^* + H.c. (13) 

The cooling laser indeed couples the vibrational and the 
electronic states of the trapped particle. 

C. Effective interaction Hamiltonian 

Let us continue by introducing an interaction picture, 
in which the Hamiltonian H in Eq. ([T]) becomes time 
independent. To do so, we choose 



(14) 



Neglecting relatively fast oscillating terms, i.e. terms 
which oscillate with 2cjl, as part of the usual rotating 
wave approximation and using the same notation as in 
Fig. [2j the interaction Hamiltonian Hi, 



Hi = ul[t,Q) {H -Ho)Uo{t,0), 



(15) 



becomes 



Hi = -hnD{ir])(T 



hg (T c+ + H.c. 
+h{A + S)a+a- +hiyb^ + hdc^c. (16) 



Here A and A-|-i5 denote the detuning of the cavity and of 
the laser with respect to the 0-1 transition of the trapped 
particle, respectively. 

In the following we assume that | A| is much larger than 
all other system parameters, 



|A| » n, \6\, g, r. 



(17) 



This condition allows us to eliminate the electronic states 
of the trapped particle adiabatically from the system dy- 
namics. Doing so and proceeding as in App.|Xj we obtain 
the effective interaction Hamiltonian 



Hi ~ hgcs D(\jf)c + H.c. -I- hv b^b + hSes c 



with (7off and Jeff defined as 
9^ 



ffcff 



2A 



and 5, 



cfi' 



A 



(18) 



(19) 



The interaction Hamiltonian Hi in Eq. ( 18 1 holds up to 



first order in 1/A. It no longer contains any atomic op- 
erators and describes instead a direct coupling between 
cavity photons and phonons. 
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D. Master equation 



These can then be used to moreover show that 



After the adiabatic ehmination of the electronic states 
of the trapped particle, the only relevant decay channel 
in the system is the leakage of photons through the cav- 
ity mirrors. To take this into account, we describe the 
cooling process in the following by the master equation 



1 



p = [Hi, p] — -K cp + pc'c) + KCpC 



4 



(20) 



with Hi as in Eq. (18), where k, denotes the spontaneous 



decay rate for a single photon inside the cavity. 



III. COOLING EQUATIONS 

In the following, we use the above master equation to 
derive linear differential equations for expectation values, 
so-called rate or cooling equations. Obtaining a closed set 
of rate equations is not straightforward due to the pres- 
ence of the displacement operator D in Eq. (10 1. To sig- 



nificantly reduce the number of rate equations which have 
to be taken into account in the following calculations, we 
first introduce two new operators x and y which replace 
the phonon and the cavity photon annihilation operators 
h and c by two new bosonic operators x and y. These 
commute with each other and provide a more natural 
description of the cavity-phonon system. 



Transformation of the Hamiltonian 



To simplify the Hamiltonian Hi in Eq. (18), we now 
proceed analogously to Ref. [SD| and define 



X — D{\ri)c. 



(21) 



This operator annihilates a cavity photon while simul- 
taneously giving a kick to the trapped particle. Since 
the displacement operator D{ir]) is a unitary operator 
(c.f. Eq. (11)) one can easily check that x fulfils the 
bosonic commutator relation 



1 



(22) 



This means, the particles created by x^ when applied to 
the vacuum are bosons. They are cavity photons whose 
creation is always accompanied by a displacement of the 



particle. Substituting Eq. (21) into Eq. (18), Hi becomes 

Hi — hgcfi X + n.c. + hScs x^ x + hi' b^'b . (23) 

In the following, we list commutator relations which can 
be derived using Eqs. (|3|, ([n]), and ([T2|, 



IT] X , 



[x, b] = ^ [x, 
t,6] = _ [x\b^ = -ifjx^ 



(24) 



\x, = —irj x{b — b^) — 
[x\b^b] = \ri{b - b^)x^ + rf x\ 



x^ X, b\ 



\x^ X. 



fe^] = [x'^x,b'^b\ =0. (25) 



Unfortunately, the operators x and b do not commute 
with each other. 

To assure that it is nevertheless possible to analyse the 
cooling process using only a relatively small number of 
cooling equations, we now introduce another operator y 
as 



y 



irj c 



(26) 



This operator annihilates phonons while simultaneously 
affecting the state of the cavity field. Using Eq. (251, 



one can show that y too obeys a bosonic commutator 
relation, 



(27) 



Using the above commutator relations, one can moreover 
show that X and y commute with each other. 



[x,y] = [x\y] 



0. 



(28) 



Notice that the above transformation of b and c in 



Eqs. (21) and (26) are unitary operator transformations 



which leave the total Hilbert space of the cavity-phonon 
system invariant. Indeed one can show that |54) 



U = exp [iryc^c (& + 6^)] 



(29) 



yields x when defining = U cW and y when defin- 

ing y a.s y = UbU^ . 

Using the x and the y operator, the interaction Hamil- 
tonian Hi in Eq. ( 23 ) can now be written as 



Hi — hgcs X + ¥[.c. + hScff x^ X + hrj'^v x^ xx^ X 

—\hrjvx^x{y~y^)+hvy^y. (30) 

This Hamiltonian is exact, since the exponential terms 
in the original Hamiltonian Hi in Eq. ( 16 ) have been 



removed via a basis transformation and not via an ap- 
proximation. 



B. Time evolution of expectation values 

In the remainder of this section, we use the interaction 
Hamiltonian Hi to obtain a closed set of cooling equa- 
tions, including one for the time evolution of the mean 
phonon number m. The time derivative of the expec- 
tation value of an arbitrary operator A, which is time- 
independent operator in the relevant interaction picture, 
equals 



(i) = Tr(Ap) 



(31) 
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When combining this result with Eq. ( 20 1 , we find that 
(A) evolves according to 



{[A,Hi]) -^niAx^fx + x^fxA) 



+K {xW{ir])AD{iriyx 



(32) 



with respect to the interaction picture which we intro- 
duced earlier in Section fll CI 

In this paper we are especially interested in the time 
evolution of the mean phonon number m which is given 
by the expectation value 



(33) 



Combining this equation with the definitions of x and y 
in Eqs. (pTl) and (p6|, we find that 



m 



^2 - ?? ki2 + rj ns . 



if we define 



7i2 = (y'^y) , ns = (x^xx^x) 
ki2 = i{x^x{y - y"^)) . 



(34) 



(35) 



This means, m and n2 are the same in zeroth order in rj. 
In order to get a closed set of cooling equations, we need 
to consider in addition the variables 



ni = {x^x) , kj = {y + y^) , 
/cio = i (y^ - y'^^) , kii = {x''x{y + y'')) 



(36) 



and 16 other expectation values which we define in 
App.[B| These are not listed here, since they appear only 
in the appendices of this paper. 



For example, applying Eq. ( 32 ) to the above introduced 



y operator expectation values, we find that their time 
derivatives are without any approximations given by 

ri2 = rjh' kii — ?]Kki2 + ri'^KHi , 

kj — 2riv ni — v kg , 

kg — v ki — 2rjHini , 

k<j = —2vkiQ + 2rjvkii+2rjnki2 — 2rj^Hini^ 

fcio = 21/ kg + 2iriv ki2 — 2rjKkii . (37) 

These five differential equations depend only on the y 
operator expectation values themselves as well as on ni, 
fell, and fci2. The time derivatives of all other relevant 
expectation values can be found in App. [C] 



C. Weak confinement regime 

Let us first have a closer look at the case, where the 
trapped particle experiences a relatively weak trapping 
potential. In this subsection, we hence assume that the 
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obtained from Eqs. (C3l, (C5|, (C7l, and (ClOl and ni, fci 



and k\2 obtained from a numerical solution of the 25 cooling 
equations which can be found in Section flll C| and in App.fC] 
Here we have r) = 0.1, v = 0.1k, S^s = 0.5 k, and g^a = 
0.1k which are typical experimental parameters for a weakly 
coupled cavity in the weak confinement regime. 



phonon frequency v is much smaller than the sponta- 
neous cavity decay rate k, while the Lamb-Dicke param- 
eter 11 is much smaller than one. 



V <^ K, and 77 ^ I . 



(38) 



When this applies, the y-operator expectation values 
evolve on a much slower time scale than all other ex- 
pectation values. In the following, we take advantage of 
this time scale separation and eliminate all relevant x and 
mixed operator expectation values adiabatically from the 
time evolution of the cavity-phonon system. The result 
of this calculation which can be found in App. [C]are ap- 
proximate solutions for rii, fcn, and fci2 up to first order 
in rj. 

Figs. [3] compares the analytical expressions for ni, fcn, 
and fci2 which we obtained in App. [C] with the results 
of a numerical solution of the full set of 25 rate equa- 
tions. For a weakly coupled optical cavity with ^off ^ 
the numerical results differ indeed only very little from 
the results in Eqs. (ICSl), (ICSl), (ICTI), and (IC10|. The ef- 



fective rate equations obtained in this subsection apply 
in this case after a short transition time of the order of 
1/k. Fig. [4] makes a similar comparison for the case of a 
relatively strongly-coupled optical cavity with g^s = n. 
In this case, there is less agreement between numerical 
and analytical results and the rate equations derived in 
this section apply really well only towards the end of 
the cooling process. Although we do not illustrate this 
here explicitly, let us mention that even less agreement 
is found when g^s ^ k. 
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FIG. 4: Difference between the analytical and the numerical 
solutions for ni, kn, and fci2 as in Fig. H] but for rj = 0.1, 
1/ = 0.1k, (Jeff = 0.5k, and g^g — k which are typical exper- 
imental parameters for a strongly coupled cavity in the weak 
confinement regime. 
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FIG. 5; Difference between the analytical and the numerical 
solutions for ni, kn, and ki2 as in Fig. [3] but for 77 = 0.1, 
V = 10 k, Sett ~ 1^, and gcs = 0.1k which are typical exper- 
imental parameters for a weakly coupled cavity in the strong 
confinement regime. 



When substituting Eqs. (Of, fC5|, (C7|, and (CIO) 



into Eq. (37), we are left with a closed set of five effective 
coohng equations which hold up to order 77^, ie. 

/ • • • • \ T T 

{h2,k7,ks,kQ,kio) ^ A/ (722, fcy, fcs, fcg, fcio) 

+ 132,^3, Pa, Psf (39) 



with 



M = 



I (2 
/ "11 






) „.(2) 



-j^ 

V oi-^l 

v(2) ^,(2) 



Ml 



^42 



^43 



V 



„(1) 



^52 



53 



44 

1v 







-2i/ 

^55 



(40) 



Each superscript indicates the scaling of the respective 
matrix element of M with respect to 77. Taking Eq. (38) 

into account, we find that the a 
good approximation given by 



,.(1) 



"12 

"42 

"52 
while 



,(i)_ 8777/52^(^2 



4<52ff) 



of M are to a very 



(k2 + 
(K2-H4<52g)2 



a 



*43 ' 



a 



a 



(2) 



33 — '^44 

327^2 Kt/ Jeff g2^ 
(«2 + 452g)2 



4y 



(2) 



4^3^ 
.,(1) 



1^14 



*55 



a 



(2) 
41 



*2+4j2^ 



647^2 7tt/ Jeff ffoff 
(k2+4<52^)2 ' 

128?/2Kt/Jeffgcff 
(K2+4j2ff)2 



(41) 



(42) 



Moreover, one can show that /3i equals, up to second 
order in 77, 



,2^„2 



/3i - 



■^v Kg. 



[(K:2+4j2ff)(7t2+4<S2ff-84fFi^) 



(«2 + 4<52ff)3 

+85e'ff (3^2 _ 4^2^)] ^ (43) 
while /32 to /^s are in first order in rj given by 



132 



k2 4- 4j2ff 



, p3 



K2+4<52ff ' 



/34 = /35 = 0. 



(44) 



We now have a closed set of five differential equations 
which can be used to analyse the time evolution of the 
y operator expectation values in the weak confinement 
regime analytically and numerically. 



D. Strong confinement regime 

In the following, we define the strong confinement 
regime as the case, where the phonon frequency v is com- 
parable or larger than the spontaneous cavity decay rate 
AC. In this subsection we therefore assume that [551 



K < v, ScS and 7/^1. 



(45) 



In this parameter regime, the time scale separation which 
we assumed in the previous subsection no longer applies. 
This means, a proper analysis of the cooling process 
should take all cooling equations into account. However, 
as we shall see below, the cooling process takes place on 
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FIG. 6: Difference between the analytical and the numerical 
solutions for ni, kn, and fci2 as in Fig. |3] but for rj = 0.1, 
1/ — 10 K, ScB — V, and g^s = k which are typical experi- 
mental parameters for a strongly coupled cavity in the strong 
confinement regime. 



a time scale which is much longer than the time scale 
given by the inverse cavity decay rate This means, 
712 evolves only on a much longer time scale than all the 
other above defined expectation values. It is therefore 
possible to simplify the 25 cooling equations introduced 
in this paper again via an adiabatic elimination. The de- 
tails of this calculation can be found in App. [D] where 
we calculate ni, fcn, and fci2 up to zeroth and first order 
in 77, respectively. 

Figs. [5] and [6] compare the obtained analytical results 
with the corresponding numerical solutions of the closed 
set 25 cooling equations. In case of a relatively weakly 
coupled optical cavity (with g^s ^ n) we find again rela- 
tively good agreement between both solutions. Although 
we now eliminate more variables from the system dy- 
namics, we find again that the results of the adiabatic 
elimination apply to a very good agreement throughout 
the whole cooling process. Less agreement is found in the 
case of a strongly coupled optical cavity with g^s = k. 
In this case, the expressions found for ni, fcn, and fci2 
apply only towards the end of the cooling process. Un- 
fortunately, it is not possible to obtain more accurate for 
this parameter regimes and the case where ^eff 3> k. This 
would require to calculate ni, fcn, and ki2 up to terms 
in the order of rf correctly which is beyond the possible 
scope of this paper. 



Substituting Eqs. (C3) and (D6) into Eq. (37), we now 



obtain only a single effective cooling equation, 



"-2 



'7c «2 



(46) 



with the constants 7c and c given by 



k2 + 4((5cff + vf 



(47) 



up to second order in 77. As we shall see below, 7c is 
the effective cooling rate for the cavity mediated cooling 
process illustrated in Fig. [T] 



IV. STABILITY ANALYSIS 



In the strong confinement regime (cf. Eq. (46 1), the 



cooling process can be described by a single effective cool- 
ing equation. Since the cooling rate 7c is always positive, 
the trapped particle always reaches its stationary state. 
However, it is not clear whether or not the same ap- 
plies in the weak confinement regime, where the cooling 
process is described by five linear differential equations 
(cf. Eq. (39 1). Proceeding as in Ref. [SO], we now have a 



closer look at the dynamics induced by these equations. 
To do so, we introduce the shifted y operator expectation 
values 



(n2, fcy, fcs, kg, fcio) 

+A./-i(/3i,/32,/33,/34,/35 



(48) 



Notice that the tilde and the non-tilde variables differ 
only by constants, namely by the stationary state solu- 
tions of the non-tilde expectation values. Substituting 



Eq. ( 48 1 into the effective cooling equations in Eq. ( 39 1 , 



they hence simplify to 



{h2,kT,ks 



kg, ho)^ 



M 772, fcy, ^8, kg, ki 



(49) 



The stationary state solution of this differential equation 
is the trivial one with all tilde variables equal to zero. In 
the following we show that the real parts of all eigenvalues 
of AI are negative, which is a necessary condition for the 
system to reach this state. 



A. Time evolution for rj = Q 



First, we calculate the eigenvalues of M in Eq. (40 1 for 
7] ~ Q and find that these are simply given by 



Ai=0, A: 



2,3 



4,5 



^2iv . 



(50) 



Taking this into account and solving Eq. ( 49 1 analyti- 
cally, we find that 



n2{t) = "2(0), 

^7(0 ^ _ ( cosvt —smut 
hit) ) \ s™^^ cosut 
kg{t) \ ( co^2vt —s\a.2vt 
koit) 



sin 2^1 cos 2^1 



kio) 

~kg{0) 

fcio(O) 



(51) 
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FIG. 7: Diagrams illustrating the time evolution of the expectation values to kio, and 722 for rj = 0.1, v = 0.1k, 
SeH = .5 K , gcff = 0.1 K and n2(0) = 100. All figures are the result of a numerical solution of the effective cooling equations 
in Eq. (39 1. In (a), only terms in zeroth order in i] are taken into account. In (b), terms in zeroth and in first order in -q 
are taken into account. In (c), all matrix elements of M in Eq. (391 are taken into account. 



These equations are illustrated in Fig. [7|a) which shows 
phase diagrams for the time evolution of the coherences 
^7 to kiQ. The fact that all points lie on a circle illustrates 
that an initially coherent state of the y particles remains 
essentially coherent throughout the cooling process. The 
numerical solution for the time evolution of 712 shows 
that, for ?7 = 0, the mean phonon number m does not 
change in time, as one would expect. There cannot be 
any cooling without an interaction between the electronic 
and the motional states of the trapped particle. 



First order corrections 



which shows a numerical solution of Eq. (49) with all 



first order corrections in taken into account. However, 
since the eigenvalues of Af have no real parts, h2 and 
therefore also the mean phonon number m, do not reach 
their stationary state solutions. Instead, h2 remains close 
to its initial value. No cooling occurs. 



C. Second order corrections 



Taking all terms in Eq. (40 1 into account, one can show 



that the eigenvalues of M are without any approxima- 
tions given by 



Calculating the eigenvalues of the matrix M in Eq. ( 40 ) 
up to first order in ij, we obtain again Eq. (50 1. All 



of them have zero real parts. But there are first order 
corrections to the eigenvectors of M. As a result, 77,2 is no 
longer constant in time. This is illustrated in Fig. [7]jb) 



Ai 

^2,3 
^4,5 



v(2) 



,.(2) 



4zy2 



,,(2)2 



„(2)„,(2) 



(52) 
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For positive effective laser detunings, the matrix element 



(2) 

a\i (cf. Eq. (42)) is always negative. This means, all 



eigenvalues of M have negative real parts, when SeS > 
0. In this case, all tilde variables are damped away on 

(2) 

the time scale given by l/aii and tend eventually to 
zero. This is illustrated in Fig. [7]jc). Now we observe an 
exponential damping of n2 which implies cooling of the 
mean number of phonons m. Analogously, one would find 
heating when solving the above equations for negative 
effective laser detunings, ie. See < 0. 

Moreover, for 6cs > 0, we find that the y coherences 
to kiQ oscillate with a slowly decreasing amplitude 
around zero. Analogously one can show that the y co- 
herences k^ to fcio oscillate with a slowly decreasing am- 
plitude around their time averages. This means, the cool- 
ing process remains stable and the trapped particle can 
be expected to reach its stationary state eventually. This 
observation is taken into account in the following section, 
where we analyse the cooling process in more detail by 
replacing the coherences fcy to fcio by their time averages. 
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FIG. 8: Logarithmic contour plot of the stationary state 
phonon number rriss in Eq. ( |56[ ) for a strongly coupled op- 
tical cavity as a function of the phonon frequency v and the 
effective laser detuning 5cff. 



V. PHONON NUMBERS AND COOLING 
RATES 

In this section, we point out that the effective cooling 
equation for 712 in Eq. ( 46 1 applies to a very good ap- 



proximation not only in the strong confinement regime 
but also in the weak confinement regime. Since n2 and 
the mean phonon number m are identical in zeroth order 
in 77, solving this equation implies that m is to a very 
good approximation given by 



m(t) = [m(0) — 
with 7c as in Eq. (|47|) and with rngg, 



7c t 



C 

7c 



(53) 



(54) 



being the stationary state phonon number for the cooling 
process illustrated in Fig. [T] in zeroth order in 77. 
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FIG. 9: Logarithmic contour plot of the effective cooling rate 
7c in Eq. (471 in units of 2g^ff/K as a function of the rela- 
tive phonon frequency i^/n and the relative effective detuning 
Scs/k for 77 = 0.01. 



The previous section shows that, in the weak confine- 
ment regime, the initial y operator coherences fcy to fcio 
oscillate relatively rapidly in time. However, since they 
oscillate with a decreasing amplitude, we can safely ap- 
proximate them by their time averages. The easiest way 
of calculating these time averages is to recognise that 
their time derivatives are equal to zero. This means, the 
time averages of kj to kio are the solutions of 



fci = for i = 7, . 



10. 



(55) 



of 772 in the strong confinement regime via an adiabatic 
elimination of fcy to kiQ. This means, the calculations in 
Section 



the wea 
tion. 



HID and therefore also Eq. (46), apply also in 



•c confinement regime to a very good approxima- 



B. Stationary state phonon number 



Exactly the same condition has been imposed in Sec- 
tion |III Pj and App.jD] when analysing the time evolution 



Substituting Eq. (47) into Eq. (54), we find that the 



stationary state phonon number 77iss is in zeroth order in 
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FIG. 10: Logarithmic plots of the time evolution of the mean phonon number m during the cooling process for typical 
experimental parameters in the strong confinement regime. The dashed lines are the result of a numerical integration of 
the 25 cooling equations which can be found in Section |III B| and in App. [C] The circles illustrate the analytical solution 
given in Eq. ( 53 1. 



r] given by 



16i^(5off 



(56) 



That this term is exactly the same as the station- 
ary state phonon number obtained by other authors 
(cf. eg. Ref. [13]), shows that our calculations are con- 
sistent with previous calculations. For example, in the 
weak confinement regime (cf. Eq. (38)), the stationary 
state phonon number TOgg assumes its minimum, when 



1 



3off 



(57) 



As already pointed out in Ref. |43j, this detuning corre- 
sponds to the stationary state phonon number 



K 



(58) 



which is in general much larger than one. In the strong 
confinement regime (cf. Eq. (45)), the stationary state 
phonon number rrigs assumes its minimum, when 



^cff 



(59) 



For spontaneous decay rates k much smaller than this 
equation simplifies to (5off — v (sideband cooling) . Substi- 
tuting this result into Eq. (56 1 and assuming k <^ we 



see that the minimum stationary state phonon number 
russ equals 

^2 

Toss = T^r^ (60) 



in this case which is indeed much smaller than one. These 
results are confirmed by Fig. [8j which shows 
function of i>/k and 6os/k. 



C. Effective cooling rate 

Let us now have a closer look at typical values of the ef- 
fective cooling rate 7c. Fig.^shows jc in units of 2g^ff/K 
as a function of i^/k and Scsji^- Since 7;. always scales as 
g^ff, the cooling rate 7c might be very small for realistic 
experimental parameters. In this case, it might seem as 
if the system reaches its stationary state, even when TOss 
is very small. 



D. Numerical results 

We conclude this section with a numerical solution of 
the full set of 25 cooling equations which we can be found 
in this paper in Section [Tll| and App.[C] Fig. 10 a) illus- 
trates the cooling process for a relatively strongly coupled 
cavity with gcff — n- Fig. [l0j[b) illustrates the cooling 
process for a weakly coupled cavity with k. We 

then compare these solutions with our analytical solu- 
tion for the time evolution of the mean phonon number 



m which takes the effective cooling rate 7c in Eq. (47) 
and the stationary state phonon number rriss in Eq. ( 56 ) 
into account. 



16j^2 



A closer look at Fig. 10 confirms that there is very 
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good agreement between analytical and numerical re- 
sults, in the case of a weakly coupled cavity. In the case 
of a strongly coupled cavity, we only observe reasonable 
agreement in the strong confinement regime when v > k. 
However, when modelling the cooling process for a weakly 
confined particle inside a strongly coupled cavity, wc find 
that the analytical expression for the stationary state 
phonon number nigs in Eq. ( 56 ) is substantially lower 



than the corresponding numerical solution. This differ- 
ence tells us that higher order terms in -q should be taken 
into account when calculating ni, fcn, and ki2 via an 
adiabatic elimination, as pointed out already in Section 



HID A much larger set of more accurate rate equations 
should be taken into account. 



VI. CONCLUSIONS 

This paper revisits a standard scenario for cavity me- 
diated laser cooling [3B1B3]- As illustrated in Fig. [T] 
we consider a particle, an atom, ion, or molecule, with 
ground state |0) and excited state |1) with an external 
trap inside an optical cavity. Moreover, we assume that 
the motion of the particle orthogonal to the cavity axis, 
ie. in the direction of the cooling laser, is either strongly 
or weakly confined and consider it quantised. The cooling 
laser establishes a direct coupling between the phonons 
and the electronic states of the trapped particle, thereby 
resulting in the continuous conversion of phonons into 
cavity photons. When these leak into the environment, 
vibrational energy is permanently lost from the system 
which implies cooling. 

As in Refs. [55H43) . we describe the time evolution 
of the experimental setup in Fig. [l] with the help of a 
quantum optical master equation. Assuming that the ex- 
cited state |1) of the trapped particle is strongly detuned 



(c.f. Eq. (17)), the system Hamiltonian can be simpli- 
fied via an adiabatic elimination of the electronic states 
of the trapped particle. We then use the resulting ef- 
fective master equation to obtain a closed set of 25 rate 
equations, ie. linear differential equations, which describe 
the time evolution of expectation values. Most of these 
expectation values are coherences. 

Since the effective cooling rate 7c (cf. Eq. ([47|) scales 



as 77 , a proper analysis of the cooling process needs to 
take terms of the order if' in the system dynamics into 
account. Instead of expanding the Hamiltonian Hi in 



Eq. ( 16 ) in 77, we solve the cooling equations for small 
Lamb-Dicke parameters rj perturbatively. The reason 
that our calculations are nevertheless relatively straight- 
forward is that we replace the phonon and the cavity 
photon annihilation operators b and c in the interaction 
Hamiltonian Hi by two new bosonic operators x and y 
(c.f. Eqs. (21) and (26)) which describe the cavity-phonon 
system in a more natural way and commute with each 
other (c.f. Eq. (28)). The operator x annihilates cavity 
photons while giving a kick to the trapped particle. The 
operator y annihilates phonons but not without affecting 



the field inside the optical cavity. 

Our results confirm that there are many similarities 
between ordinary and cavity mediated laser cooling |43| . 
However, for a weakly confined particle inside a strongly 
coupled cavity, a comparison between analytical and nu- 
merical results suggests that more detailed calculations 
are needed to model the cooling process accurately. Our 
analytical calculations are designed such that they cal- 
culate Toss in zeroth order in the Lamb-Dicke parameter 



•q (cf. Eq. (56)). This means, we neglect higher order 
terms in -q in the rate equations, whenever possible. Our 
numerical calculations however take all available terms 
in the above mentioned 25 rate equations into account. 
The difference between analytical and numerical results 
means that terms of higher order in rj are not negligible, 
although this might seem to be the case. Unfortunately, 
calculating mgs systematically up to first order in ry, ei- 
ther analytically or numerically, would require to take 
considerably more than only 25 cooling equations into 
account. 

The above observation is nevertheless interesting, 
since the cooling of a weakly confined particle inside a 
strongly coupled cavity is of practical interest for the 
cooling of molecules. Realising a very strong coupling 
between a trapped particle and the field inside an optical 
cavity is in principle feasible [131 [2]. Over the last years, 
experiments have been performed with a continuously 
increasing ratio between the cavity coupling constant g 
and the spontaneous cavity decay rate k. Even larger 
ratios g/n are expected to occur when trapping large 
molecules, whose electric dipole moment Dpi can be 
much larger than that of an atom, inside an optical 
cavity. Such molecules can experience a relatively large 
effective cavity coupling constant t^cff- 
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Appendix A: Adiabatic elimination of the electronic 
states 



In the following, we write the state vector of the atom- 
cavity-phonon system as 



1 00 00 



J— m— n— 



(Al) 



where |j) and |to) denote the electronic and the vi- 
brational energy eigenstates of the particle and where 
|n) is a cavity photon number state. According to the 
Schrodinger equation, the time evolution of the coefh- 
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cient Cj'rn'n' is given by 
1 

1 



J. 'JO LXJ 

- ^ ^ ^ Cjmn {j'm'n'\Hi\jmn) 



j—0 m—Q n—0 

(A2) 

Given condition (171, the coefficients Cj'm'n' with j' = 1 
evolve on a much taster time scale than the coefficients 
with / = 0. Setting the time derivatives of these coeffi- 
cients equal to zero, we find that 



Clr. 



1 



2A ^ ^ 

m—O n—0 



COn 



x{m'n'\ (nD^i-q) +2gc) \mn) , (A3) 

if the particle is initially in its ground state. This equa- 
tion holds up to first order in 1/A. Substituting this re- 
sult into Eq. ( A2 ) and neglecting an overall level shift, we 



obtain the effective interaction Hamiltonian in Eq. ( 18 ) 



Appendix B: Relevant expectation values 

The calculations in Apps. [C| and [P] require in addition 



to the expectation values defined in Section IIIB the x 
operator expectation values 



ki = {x + x') , k2 = i{x 

fcg EE (x^ + x'' 

k^ = {x\x A- x^)x) , fcg = i [x^ [x — x^)x) . (Bl) 

Moreover we employ in the following the mixed operator 
expectation values fci5 to k22 which are defined as 



fcl3 


^ 




fci4 


^ i((^ 




fcl5 




x'^){y~y'')) , 


fcl6 


^ i((^ 


+ x^){y-y^)), 


fci7 




x^){y + y^)}, 


fcl8 


^ i((^ 


^x^)iy + y^)), 


fcl9 


^ ((^- 


x^){y'-y^')), 


^20 


^ i((x 


+ x^){y^-y^')), 


^21 




x^){y' + y^')), 


^22 




-x^){y'+y^')). 



(B2) 

Their time derivatives of these and other expectation val- 
ues can be found in App.[C| where they are use to obtain 
a reduced set of effective cooling equations. 



Appendix C: ni, fcn, and fci2 in the weak 
confinement regime 

In this appendix, we derive approximate solutions for 
the expectation values kn, ki2, and ni for the weak con- 
finement regime (cf. Eq. (38)). This is done via an adia- 
batic elimination of the x and the mixed operator expec- 
tation values which all evolve on the relatively fast time 



scale given by the spontaneous cavity decay rate n. To 
indicate the scaling of variables, we adopt the notation 



X = + X 



(1) 



r(2) 



(CI) 



The superscripts indicate the scaling of the respective 
terms with respect to rj. As we shall see below, the ex- 
pectation values fell, fci2, and ni need to be calculated 
up to first order in rj. Let us first have a look at 



^^2^ ^-nd rif\ 



Using Eq. (32 1 and setting 77 = 0, wc find that rii, n^^ 



and fci to fcg evolve in zeroth order in 77 according to 



ni 

ki 

k 

k 
ki 



5cff k2-K.n1, 

5off (^2 + 2/05) + n {ni~ 2^3) , 

2.gcff + (5cff ki - -Kk2, 

-Sgcff - 2(5eff ki - nkz, 
2gcff ki -f 2(5eff ks- nki, 

3 i. 
^i^k5, 



5off ki - Scff ke 
gcs (4rii - ks) 



(C2) 



These equations form a closed set of differential equa- 
tions. Eliminating the above x-operator expectation val- 
ues adiabatically from the system dynamics, we find for 
example that ni is in zeroth order in 77 given by 



.(0) - 



eff 



k2 4- 4(5,2 



(C3) 



off 



In addition we obtain expressions for kf^ , ki^^ , k^^ , 



and n 



(0) 



2 I "-5 ' ""6 ' 

These are used later on in this appendix to 



calculate k[^^ and k[^2 ■ 

Setting r/ = and using again Eq. ( 32 1 , we moreover 
find that the time evolution of the mixed operator coher- 
ences kii and ki2 and fcis to /cig is in zeroth order in rj 
is given by 



fell 
ki2 

ki5 
he 
kn 



5cff ^18 - i^ki2- nkii , 
-9cs ki5 + vkii - Kki2, 

-2.gcff ^8 - Scs kie - v fcig 

(Joff fcl5 -I- i^fcl7 - ^Kfcie , 



-K,k 



15 : 



-(5cff /cis -vki 



2gcff k^ + (5cff fcir + v fci5 



-K fci 



(C4) 



These six equations too form a closed set of cooling equa- 
tions which describe a time evolution on the time scale 
of the spontaneous cavity decay rate k. Taking this into 
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account, eliminating fcn and ki2 and ^15 to fcig adiabat- 
ically, and neglecting terms proportional to which are 
much smaller than the remaining terms, we find that 



'^11 — 



k7- 



^2 + 4(52fj 



(C5) 



In addition we obtain expressions for kf^ and k^^ which 
are used below in the next paragraph. 

Proceeding as above but taking terms up to first order 
in ry into account we find that the first order contributions 
of the X operator expectation values ni, ki, and /c2 in 
Eq. ( 35 ) evolve according to 



^(1) 



k['^ 



h.2 



5ofr ^2^' 



= -Scffk. 



= (Scff k[^^ - rjiy k'^^g' - ^nk. 



,(1) 



(0) 



1 



(C6) 



These equations form a closed set of cooling equations, 
when the above mentioned results for k^^ and k^^ are 
taken into account. Eliminating ni, ki and ^2 adiabati- 
cally and neglecting all terms proportional to , wc find 
that 



,(1) 



ks. 



(C7) 



This means, n'^^ follows the time evolution of fcg adia- 
batically. 



In order to calculate k[^^ and k\2 , we need a closed 
set of cooling equations which applies up to first order in 
r] correctly. Applying Eq. (32 1 again to fcn and ki2 and 
ki5 to fcig, we find that the time derivatives of their first 
order corrections in rj are given by 



(1) 



'^11 ~ i/cff "-IS 



(1) 



'^12 



'^16 



'^17 



-5cff k[l^ - 
+2r]K 
-2riKkf^ 



"-12 



(1) 

18 



- 2r]K 



(0) 

(0) , 



(0) 



2k 



13 



- K k 

JO) 
^21 



(1) 



2'*'''15 ' 



1 



2^*^16 



(1) 



JO) 



- 2A;^4'' fc22 



(0) 



--Kk^'^ 



k 



(1) 

16 



JO) 



2k 



(0) j,(0) 



<5cff' k\y + V k^^^ + r\v 



JO) 



2/c, 



(0) 1.(0) 



6 



20 



-K fc 



(1) 
18 



(C8) 



Substituting the definitions of the mixed-particle expec- 
tation values /cia and fci4 and fcig to fc22 into Eq. ( 32 ) and 
setting 77 = 0, we moreover find that 



^13 = 

fci4 = 

^19 = 

^20 = 

ki\ = 

ki2 = 



-Setfku - ^Kfcl3 , 

2gcS n2 + Scff fcis - -nku, 

-2gcS kio - Sctf ^20 - 2t/ fc22 - -k fcl9 : 
<5off ^19 + 2zy fc21 - ^20 , 



-Sctf k22 - 2u k20 - -K ^21 , 



1 



2gcS kg + Scs k2i + 2u kiQ - ^22 ■ (C9) 

These final six differential equations hold in zcroth order 
in 77. Setting the right hand side of these and of the 
cooling equations in Eq. (C8) equal to zero, we finally 
obtain the expressions 



(1) _ 



"■11 



Ji) 

"12 



f 2 , AX2 \2 ^'"off '^lO + H 
^«(«2+4<52g)4 
(k2 -t- 45^ff )2 

32,7g4^(3«2-4<S2g) 



Again we neglected terms proportional to ly^ , since these 
are in general much smaller than the remaining terms. 

Appendix D: ni, fcn, and ki2 in the strong 
confinement regime 

Let us now have a closer look at the strong confinement 
regime (cf. Eq. (45)). However, different from the pre- 



vious subsection, we no longer assume that some system 
parameters are much smaller than others. The reason 
that we nevertheless obtain relatively simple expressions 
for the quasi-stationary state solutions for ni, fcn, and 
ki2 is that we eliminate in the following not only the a; 
and the mixed operator expectation values, but also the 
y operator coherences ky to kiQ. From Eq. (37) we see 



that calculating 77.2 up to second order in 77 requires know- 
ing ni in zeroth order in 77. Having a closer look at the 
above cooling equations, we see that the expression for 
in the strong confinement regime is the same as the 



JO) 



expression in Eq. (C3). In addition, wc need to calculate 



fell and ki2 up to first order in 77. 

Using again Eq. (37), setting 77 = and eliminating 



the y operator coherences adiabatically from the system 
dynamics, we find that kj to kio all equal zero in zeroth 
order in 77, 



k 



(0) 



JO) 



JO) 



JO) 
"10 



0. 



(Dl) 
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Taking this into account when ehminating the mixed 
operator expectation values fcn, ki2, and /ci5 to fcis in 



Eq. ( C4 ) adiabatically, we now find that all of them van- 



ish in zeroth order in 77, 
'^11 - 



^12 



0. 



(D2) 



This means, the time derivative of n2 in Eq. (37) scales 
as r]^, at least to a very good approximation. 

To calculate ku and ki2 up to first order in ij, we have 
agai n a c loser look at Eq. (37). Using this equation and 
Eq. (C3), one can show that the y-coherences k^ and fcg 



are in first order in r/ given by 



.(1) 



cff 



.(1) 



:2 + 4,52 



(D3) 



cff 



Using Eqs. (C9| and (Dl), we see in addition that 

4K.goff 



, (0) _ 8^cff.9cff 



u(o) 

^19 



k 



(0) 
20 



4(^2 

l(0) _ 
'^21 - 



l(0) _ 
^14 — 
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while k[\^ , k[^2 j ^le'' : ^-nd fcjy^ evolve as stated in Eq. (|C8 1 
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Substituting Eqs. (D3| and (D4) into these equations, 
using the solutions for n^*^-* , ng^and the coherences k\ 
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and 



the system dynamics, we obtain 
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adiabatically from 
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Applying Eq. (32) again to ku, ki2, and fci5 to fcig, we 



find that the time derivatives of the fcis to kis in first 
order corrections in rj are now given by 



with the constant ^"^ defined as 
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